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ABSTRACT 

We revisit the metal budget at z ~ 2, and include the contribution of the intergalactic 
medium. Past estimates of the metal budget indicated that, at redshift z ~ 2.5, 90% 
of the expected metals were missing. In the first two papers of this series, we already 
showed that ~30% of the metals are observed in all z ~ 2.5 galaxies detected in 
current surveys. This fraction could increase to < 60% if one extrapolates the faint 
end of the LF, leaving >40% of the metals missing. Here, we extend our analysis to the 
metals outside galaxies, i.e. in intergalactic medium (IGM), using (1) observational 
data and (2) analytical calculations. Our results for the two are strikingly similar: 
(1) Observationally, we find that, besides the small (5%) contribution of DLAs, the 
forest and sub-DLAs contribute subtantially to make <30-45% of the metal budget, 
but neither of these appear to be sufficient to close the metal budget. The forest 
accounts for 15-30% depending on the UV background, and sub-DLAs for >2% to 
<17% depending on the ionization fraction. Combining the metals in galaxies and in 
the IGM, it appears now that >65% of the metals have been accounted for, and the 
'missing metals' problem is substantially eased. (2) We perform analytical calculations 
based on the effective yield-mass (y e s-V c ) relation, whose deficit for small galaxies 
is considered as evidence for supernova driven outflows. As a test of the method, 
we show that, at z = 0, the calculation self-consistently predicts the total amount 
of metals expelled from galaxies. At z — 2, we find that the method predicts that 
25-50% of the metals have been ejected from galaxies into the IGM, consistent with 
the observations (<35%). The metal ejection is predominantly by Lb < \L* B (z = 2) 
galaxies, which are responsible for 90% the metal enrichment, while the 50 percentile is 
at L ~ jqL* b (z = 2). As a consequence, if indeed 50% of the metals have been ejected 
from galaxies, 3-5 bursts of star formation are required per galaxy prior to z — 2. The 
ratio between the mass of metals outside galaxies to those in stars has changed from 
z = 2 to z = 0: it was 2:1 or 1:1 and is now 1:8 or 1:9. This evolution implies that a 
significant fraction of the IGM metals will cool and fall back into galaxies. 
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1 INTRODUCTION 

Locally, our picture of the metal (and baryon) budget has 
become more and more complete over the past few years 
i|Fukugita et all 1 19981 : IFukugita fc Peebles! l2004f ). Roughly 
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speaking, 30 % of the baryons are in the Lya forest 

ilStocke et al.ll2004j l. 50% are in a warm-hot phase (WHIM, 
iTripp et all 12004 ), 5-10% are in the intra-cluster medium 
(ICM), and 10% are in stars. So, even-th ough WHIM (e.g. 
ISembach et alj|2004 iNicastro et "ail 120051 ) and intra-group 
medium I Fukugita fc Peeblesl 2004 ) dominate the baryon 
budget, they contain a minor fraction (< 10 %) of the met- 
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als. Indeed, about 80-90% of all the metals produced by 
type II and type la supernovae are locked in stars (10%) 
and stellar remnants (80%) such as white dwa rfs, neutron 
stars and black holes (IFukugita fe Peeb les 2004). 

Ten Gyrs ago (z ~ 2.5), the situation was very dif- 
fer ent;_^lost_£90%)_oft h^_baxyons_were_J^i the Lya forest 



Rauch et all Il997l: IPenton et all l200d : ISchavd l2001al : 



ISimcoe et aT1l2004h . but our knowledge of metal abundances 
is still highly incomplete. In the past, it was realized that 
only a small fraction (20%) of the expected metals is ac- 
counted for when one adds the contribution of the Lya for- 
est (-Wjjj = 10 13-17 cm -2 ), damped Lya absorbers (DLAs) 
(iVjj : > 10 20 ' 3 c m" 2 ), and galaxies s uch as Lyman break 
galax ies (LBGs) l|Pettini et all Il999l ; IPaeell 120021 ; IPettinil 
l2003fl . 

As discussed in iPettini et ail (|2003h . either the 
missing metals were expelled into the intergalactic 
me dium (IGM) via galactic w inds (as already discussed 
by"" 



Larson fe Dinersteinl Il975 | ) since LBG s drive winds 



(AaciDcrger ex ai. ziuuj; anapiev et ai. 
one seen locally dLehnert & Heckman 


zuuj j mucn iikc tnc 
1996a; Dahlem et al. 


1997 


; iLehnert et all 19991; Heckman et al. 


20001: Martin 



1999), or they are in a galaxy population not accounted for 
so far. 

In lBouche et all l|2005l ) (hereafter paper I), we showed 
that only 5% (and < 9%) of the ex pected metals are in 
submm selected galaxies (SMGs) and in lBouche et all (|2006l ) 
(hereafter paper II), we revisited the missing metal problem 
in light of the several new 2 ~ 2 galaxy populations, such 
as the 'distant red galaxies' (DRGs) and the 'BX' galaxies. 
Paper II showed that the contribution of z — 2 galaxies 
amounts to 18% for the star forming galaxies (including 8% 
for the rarer AT-bright galaxies with K s < 20) and to 5% for 
the DRGs (J — K > 2.3). Thus, adding the contribution of 
star-forming 'BX' galaxies, DRGs, and the SMGs, the total 
contribution of 2 ~ 2 galaxies is at least 18 + 5 + 5 ~ 30% of 
the metal budget. As shown in paper II, if one extrapolate 
these results to the faint end of the luminosity function, 
<60% of the metals are in 2 ~ 2 galaxies, leaving > 40% 
missing. 



The goal of this paper is to use observational data (sec- 
tionf2| and analytical calculations (sectionfSJ) to gain insights 
on the missing metals. In section [2] we include the contribu- 
tion of metals in various gas phases, traced by QSO absorp- 
tion lines. The dominant contributors are the forest (with 
logiV Hl = 10 13 ~ 17 ) and the sub-DLAs with logiV Hl = 
10 19-20 ' 3 . We will show that neither the forest alone, nor sub- 
DLAs, can account for the remaining 40% of missing met- 
als. In section [3] we explore the possibility that much of the 
remaining missing metals were ejected from galaxies based 
on 'effective yield' arguments. Namely, locally small galax- 
ies have lost a significant fraction of their metals, reflected 
by their offset from the closed-box expected yield (e.g. 
Pilvugin fe Ferrinilll998l ; lKoppen fe Edmunds! 19991 ; iGarnettl 



20021 ; lDalcantonll2006l l Indeed, if all galaxies were evolving 
as 'closed boxes', their metallicity Z would be inversely pro- 
portional to the gas fraction p, i.e. Z cx — \n(p), and 
they would all reach approximatel y solar metallicity once 
they convert their gas into stars (e.g . Scarlc & Sargcn fc[l972t 
Audouze fe Tinslevj[l976l; iTinslev feLarson 1978; Ed munds! 



tionality constant, is defined as; 
Z 



VcB = 



In(l/A*) 



(1) 



The ratio between the solar yield and y e g gives the fraction 
of metals that were ejected (see section 3.1). 

In the remainder of this paper, we used a Ho = 
70 h 70 kms -1 Mpc~\ Q M = 0.3 and fi A = 0.7. We use 
a Z© = 0.0189 and a Salpeter IMF throughout. 



2 THE METAL BUDGET 

2.1 The Metal Production Rate and Assumptions 

Stellar nucleosynthesis and star formation govern the pro- 
duction of metals. Some of the metals remain locked up 
in stellar remnants or long-lived stars, some are returned 
into the interstellar medium (ISM) and another fraction 
is expelled into the intergalactic medium (IGM) via galac- 
tic winds and outflows. Under the assumption of instanta - 
neous recycling approximation Q (Se arle fe Sargent] I1972T ). 
the global metal production rate pz can be directly re- 
lated to the star formation rate den sity (SFRD) p+ through 
the IMF-weighted yield < p. > (ISearle fc Sargentl Il972l ; 
ISongaila et allll990r , livladau et al1ll996h : 



pz = <Pz > p* , with 
< p z > = / dm p z (m)m<j)(m) 



(2) 
(3) 



where <j>(m) is the IMF, and p z (m) is the metal yield for stars 
of mass m. The total expected amount of metals pz, eX pected 
formed by a given time t is the integral of Eq. [2] over time. 

T he IMF-weighted y ield in Eq. [2] is equal to or 
2.4% (|Madau et all Il996l ) using a Salpeter IMF and the 
type II stellar yie l ds (fo r solar metallicity) p z (m) from 
IWooslev fe Weaverl (|l995l )Pl It should be noted that chang- 
ing the IMF (or the low mass end) will change both p\ and 
p z in the same way leaving pz unchanged since p\ is mea- 
sured from a stellar luminosity density pl- In other words, 
the metal production rate is directly related to the mean 
luminosity density pL via pl = e pz, where e is indep endent 
of the IMF jSongaila et al.lll990l ; iMadau et al1ll996l) . 

The SFRD in Eq.f2] p*, is inferred from rest-frame UV 
surveys and at redshifts greater than 1, appears to be con- 
stant from a redshift of 2 ~ 4 to z ~ 1 at a level of about 
p\ ~ 0.1 /170 Mq yr -1 Mpc~ 3 once corrected for dust extinc- 
tion (iLillv et al|l996l;lMadau et al.ll996l ; |Steidel et al.lll999l: 
Dickinson et all I2OO3I; iGiavalisco et all [20041; iDrorv et"al 



Beacom 2006; Pantcr et al 



1990). The effective yield, y c a, which is simply the propor- 



2005; Hopkins 2004; Hopkins , 
l2006l ; lFardal et alj|2006h . 

In order to compute pz(t) from Eq.[2] we parameterize 
the SFRD, p\, as shown in Fi g. [T] (bottom). Th e solid line 
shows p\ (t) parameterized by ICole et ail (|200ll ) for an ex- 
tinction of E(B - V) = 0.10. The dashed line shows p\(t) 



1 The instantaneous recycling approximation applies here since 
the timescale for massive stellar evolution (10 7 yr for a M ^ 
15 Mq star) is much shorter than cosmic timescales (10 9 yr) at 
z ~ 2.5. 

2 The c ontribution from ste llar winds from massive stars is neg- 
ligeable l lHirschi et al.ll2005h . 
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set to 0.1 /170 Mq yr 1 Mpc 3 above redshift 2 = 1, and lin- 
early proportional to z below z = 1. The dotted line shows 
p*(t) set 0.15 h 7 o M Q yr _1 Mpc" 3 at high redshifts, while 
keeping the same decline below z — 1. 

We first check that the integrated SFRD gives the 
amount of stars using p* = (1 — R) ■ f Q dt p*(£) , R is 
the mass fraction r ecycled into the ISM. For a Salpeter 
IMF, R equals 0.28 (|Cole et al.ll200lh . The top panel shows 
p*(t) for the three parameterizations. Both the solid and 
dashed curves reach t he local stellar m ass density shown by 
the sh aded area for ([Cole et a.1.1 l200lh and the filled sym- 
bol for iBell et al.l (l2003h l 3 l. The agreement was also found 



0.2 0.5 



by others (e.g. 



iRudnick et al.,1 



Rudnick et aLll2003l ; iDickinson et~ai] 120031; 

2006). Recent studies pointed to some ten- 



sion between the local p* ( z = 0) and the int egrated SFRD 
jHopkins fc BeacorrJliooll iFardal et allbOOfJ ). 

The overlap between different galaxy populations is 
very much unknown. With respect to the SFRD, two 
populations, LBGs and SMGs, reaches sim ilar levels (e.g. 
ISteidel et al] Il999l ; IChapman et"aH 120051 ). after a dust- 
correction for LBGs of a factor of ~ 5. Should the SFR 
of these two populations be added, which would double the 
SFRD at z = 2? The dotted line in Fig. [Qsliows that, if one 
doubles the SFRD at z = 2, one over-predicts the observed 
stellar density p*(z = 0) by ~ 40% if the SFRD is increased 
by 50%. Thus, there is little room to add the SMG contribu- 
tion to the SFRD at high-redshifts (~ 0.1 M Q yr -1 Mpc" 3 ), 
given the large dust-correction applied to LBGs that al- 
ready encompasses various galaxy populations. This con- 
clusion is very general and does not depend on the tension 
be tween p+(z = 0) and the in t egral of the SFRD discussed 
in Irlopkins fc Beacoml i|2006r ); IFardal et all (|2006r ), adding 
the SMGs, i.e. doubling the SFR at z ~ 2, would increase 
the discrepancy further. 

2.2 The expected amount of metals 

The right axis of FigQJtop) shows the metal density pz{i) 
found from integrating Eq.[2]over time. From this figure, one 
sees that at z = 0, the amount of metals formed (by type II 
supernovae) is expected to beQ 



/OZ,cxp.(z = 0) 



2.13 x 10 7 M Mpc" 3 , (4) 



if we integrate the SFH from t = 13.7 Gyr until the present. 
As a cross -check, this number ( Eq.ffl is consistent with the 



findings of lFukugita fc Peebles! (120041 ). Indeed, their Table 3 
list the amount of metals in different classes of objects. Given 
that the mean yield < p z >= 1/42, used in Eq. 3] is calcu- 
lated including only stars with m > 10 Mq, we exclude 
the contribution from main se quence (MS) stars and from 
white dwarfs from the table of lFukugita fc Peebles! (|2004h . 



3 We converted the results of lBell et all J2003T) to a Salpeter IMF 
by multiplying their results by xl.6. 

4 In the literature, it is so metimes useful to e xpress pz relative 
to the baryon density (e.g. |Pettinlll2003l , lioOtih : ~Z b = p z /p b /Z Q . 
This quantity represents the 'fraction of baryons with solar metal- 
licity' or the 'mean metallicity' of the universe if the metals 
were spread uniformly over all the baryons. Eq. [4] corresponds 



to Z h 



■ 0.20 



<Pz> 0.0189 
1/42 



, which is close to the mean metal- 
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Figure 1. Bottom panel: Star formation history (S FH) of the 
univer se. The solid line shows p*(t) parameterized bv lCole et al.l 
ll200ll) for a extinction of E(B - V) = 0. 10. The dashed line shows 
p*(t) set to 0.1 /170 Mgyr -1 Mpc -3 beyond redshift z = 1, and 
linearly proportional to z below z = 1. The dotted line shows p* (t) 
set 0.15 /170 Mq yr -1 Mpc -3 at high redshifts, while keeping the 
same behavior at low redshifts. Top panel: Mass stellar build-up 
(p*(i)) for the same three parameterizations. The left axis shows 
the visible stellar density ((1 — R) ■ J SFH), where R is the re- 
cycled fraction, while the right axis shows the metal production 
(< p z > • J SFH). Both the solid and dashed curves reach the 
observational constraint of p+(z = 0) ( Cole et al. 200l|), shown by 
the shaded area. The dotted line clearly overshoots the observa- 
tions, leaving little room on the SFR at high-redshifts . As a conse- 
quen ce, the contr ibutions of LBGs (dus t-corrected, ISteidel et al.l 
1 19991) and SMGs dChapman et al.ll2005T) cannot simply be added. 



and find that the z = metal density is pz, bscrvod (2 = 0) ~ 
2.3 x 10 7 M Q Mpc" 3 , which is quite close to Eq.g] In the re- 
mainder of this paper, we will refer to these numbers as 'the 
z = metal density', but one should keep in mind that it 
does not include the contribution from white dwarfs (SN la) 
and locked in the MS (see Table [Cl|l since we will focus on 
the type II yields at z ~ 2.5 throughout this paper. 

At high-redshifts, if we integrate the SFH, from z = 4 
to z = 2, corresponding to a time interval of 1.68 h^o Gyrs, 
the expected amount of metals formed is 



pZ,cxp.( z = 2) 



4.0 x 10 



0.035- 



6 P* <Pz 



0. 1 1/42 

1, > 0.0189 



> h° 70 M 



Mpc" 3 (6^ 



1/42 



-1 5 



By comparing Eq. [5] with Eq. [4] and from Fig. [T] one 
sees that about 1/4 or 1/5 of the metals are alread y pro- 
duced by z = 2 (e.g. iPagell |2002| ; |Pettini|[200l . l2006h . This 
is consistent with the stellar mass density evolution stud- 



licity of the ICM (Rcnzini 2004, and references therein). 



5 We note that iPettinil ll2006h uses p z , 
10 6 M Mpc" 3 < p z >= 1/64, and Z e 
Z h = 0.045. 



cxp. — 3.4 x 
0.0126, which gives 
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ies (e.g. iRudnick et al.J|2003l , l200d ; |Panter et al.ll2006l) that 
have found similar amount of evolution. 

If we integrate the constant SFH (at 
0.1 M yr _1 Mpc" 3 ) from z = 10 to z = 2, then the 
amount of metals formed by z — 2 is 

,6 < Pz > . 



pz,cx P .( z > 2) 



6.4 x 10° 



1/42 



(6) 



2.3 Metals in galaxies 



Some argued that the dusty ISM of submm galaxies could 
harbor most of the remaining missing metals. However, in 
paper I, we showed that the remaining missing metals can- 
not be in SMGs based on direct metallicity and gas mass 
measurements. Only 5% (and < 9%) of the expected metals 
are in SMGs. 

The second paper of this series (Paper II) showed that 
the contribution of z = 2 galaxies amounts to 10% for the 
star forming 'BX' galaxies, to 8% for the rarer A'-bright 
'BX' galaxies with K 3 < 20, and to an additional 5% for the 
DRGs (J-K> 2.3). 

We note that the numbers quoted in paper II for the 
BX galaxies were technically based on their stellar mass es- 
timates (using a Salpeter IMF) combined with their metal- 
licity, i.e. without the contribution of the ISM which was 
unconstrai ned at the time. However, based on Ha flux mea- 
surement, lErb et all (|2006l ) estimated the gas fractions of 
'BX' galaxies to be ~ 50%. Thus, the contribution of the 
'BX' galaxies alone may be as high as 15%. DRGs are gener- 
ally likely to be poorer in gas. But, given that stellar masses 

of high-redshift galaxies are known to be overestimated 
(for a Salpeter IMF) when one compares to their dynam- 
ical masses Md yn (see lForster Schreiber et al.j|2006h . our es- 
timates in paper II can be viewed as inclusive of all the 
baryons (gas and stars). 

All in all, the total contri bution of the known z ~ 2 
populations is ~ 30% (see also iPettinil l2006h of the metal 
budget corresponding to a cosmic metal density of 



PZ.gala 



1.1 x 10 6 M Mpc' 



(7) 



Furthermore, since most high-redshift surveys detect 
galaxies to a luminosity comparable to L* , in paper II, we 
estimated the contribution of the fainter galaxy population 
using a metallicity-luminosity relation. We found that the 
< L* galaxy population could double this sum. In other 
words, currently known galaxy populations at z — 2 can only 
account for 30 to 60% of the metals expected, and > 40% 
are unaccounted for (see summary in Table [TJ. 



2.4 Metals in absorption lines 

We separate our analysis of the metals in the IGM accord- 
ing to the H I column density. The forest with iVjj j < 
is discussed in section 12.4.11 sub-DLAs with 
10 ly < iVur < 10 20 ' 3 cm" 2 in section I2A21 and DLAs with 
iV Hl > 10 4d ' : 



> L9 < Nkt < 10* 

! ' 3 cm -2 in section [231 



2.4-1 Metals in the forest 

In the low-density IGM (with JV Hl < 10 17 cm" 2 ), C, Si, 
and O, as traced by Civ, Siivand Ovi, are useful probes 



of the metals contained in the IGM ( e.g. Son gaila 



Aguirr e et al 



Simcoe et al. 



.20021; ISchave et alJ |2003 ; iPettini et al J 120031 : 
20041 : lAguirre et alj|2004 ISongailall2()05h 



2001 



Several approaches have used to estimate the metal con- 
tent of the IGM using these ions. For instance, the ratio 
of C iv to H i usin g the pixel optical depth method (e.g. 
lAguirre et all [2002) can be converted into carbon abun- 
dances given an ionization correction (e.g. ISchave et al.l 
2003). These can be computed as a function of density n-pj j 
for a given UV background (UVB) model using codes such as 
CLOUDY. Hydrodynamical simulations can provide inter- 
polation tables of the density (and temperature) as a func- 
tio n of the Lyq o ptical depth (and redshift). This was done 
by ISchave et al.l (|2003l ). under several assumptions regard- 
ing the UVB model use d in generating th e ioniz ation cor- 
rections. For a UVB from lHaardt fc Madaul (120011 ) including 
quasars and galaxies (hereafter 'soft UVB'), they found 



He 



2.3xl0- 7 10[ C / H ]+ 2 ' 8 (J±-) 
V 0.0457 

3.1 x 10 4 1 [ C / H ]+ 28 



pctfGM ^ 3.1 x 10- I0 L ~'" J — ' (8) 

This corresponds to a metallicity contribution of Pz.igm — 
3 x 10 5 M Q Mpc" 3 , or only about 8% of the metal budget 
This study had a threshold at -rjjj ~ 1, corresponding t o 
7V Hl ~ 10 135 cm" 2 . The work of ISimcoe etail [|2004h . 
based on line fitting and on different assumptions, gave sim- 
ilar results ([C/H] =— 2.2) and would be almost identical 
([C/H] =— 3.1) under similar assumptions about the UVB. 
This study had a threshold at iVrj t > 10 14 cm" 2 . 

As f or silicon abundances , the firmest estimates are pro- 
vided bv lAguirre et all i|2004h , who studied the forest metal- 
licity by analyzing Si IV and C IV pixel optical depth derived 
from high quality Keck and VLT spectra. They find that 
[Si/C] ranges from [Si/C] ~ 0.25 (for a very soft UVB) to 
[Si/C] ~ 1.5 (for the 'hard' UVB). For a fiducial 'soft' UVB 
model, they fit a value of [Si/C] = 0.77±0.05 for gas in over- 
densities of 5 > 3 or A^jjj > 10 14 cm" 2 . This then gives 
a metallicity contribution of [Si/H] 
sponding to 



-2.03 ± 0.14, corre- 



Si, IGM 



3.4 x 10" 7 l [ Si / H ]+ 2 



0.045 



psi.iGM ^ 4.6 x 10 4 lo[ Si / H l +2 - M Mpc" 3 . (9) 

Sin ce about 3-5% of type II supernova metal production is 
Si (|Samlandlll998h , this corresponds to 



Pz.igm =i 1.16 x 10 6 ^ Mq Mpc" 3 , 
~Z b ~ 0.010, 



(10) 



i.e. ~30% of the metal budget. 

This result is highly sensitive to the hardness of the 
assumed UVB, since softening the UVB (for example) lowers 
both the inferred [C/H] and the inferred [Si/C]. However, the 
UVB hardness has the opposite effect on [O/C] and [O/H], 
so they are also quite usefu l to examine. 

Using a 'hard' UVB, ISimcoe et al.l (|2004h measured 
[C,0/H] from which they infer flo ^ 1.4 x 10" 6 h' 1 ~ 
2.0 x 10" 6 using [0/C]~ set by the UVB. The corre- 
sponding density of Oxygen is po ~ 3.4 x 10 5 M Mpc" 3 , 



p z , forost ~ 4.6 x 10 5 -^- M Mpc" 
0.60 



(11) 



The missing metals problem. Ill 5 



Table 1. The metal budget for z ~ 2 galaxies. 



Class 


Pz 


n z 


z b 


Pz/PZ.tot a 


Ref b 


Note 






(Mq Mpc" 3 ) 




(z e ) 


(%) 








SMGs 


< 3.6 x 10 5 


< 2.7 x 10~ 6 


<0.0031 


<9 


1 


> 3mJy (Z = 


2Z Q ) 


SMGs 


1.8 x 10 5 


1.6 x 10" 6 


0.0015 


5 


1 


> 3mJy (Z = 


\Zq) 


BX 


3.8 x 10 5 


2.8 x 10~ 6 


0.0033 


10 


2 


> L* 




BX+K <20 


3.0 x 10 5 


2.2 x 10~ 6 


0.0027 


8 


2 


> L* 




DRGs 


2.0 x 10 5 


1.5 x 10" 6 


0.0018 


5 


2 


J - K > 2.3 




Total Observed 


1.1 x 10 6 


7.8 x 10~ 6 


0.010 


~30 


2 


> L* 




Total Inferred 


2.2 x 10 6 


1.6 x 10" 5 


0.020 


~60 


2 


all L 




Missing 








> 40 




all L 





"The fractional contributions are calculated using the amount of metals expected from the SFH: pz,tot = 4 X 10 6 Mq Mpc 3 (Eq.[5|. 

'References: (1) Paper I, (2) Paper II. 



Table 2. Metals in the forest. 



Element 


Px 

(Mq Mpc" 3 ) 


n x 




Pz 

(Mq Mpc" 3 ) 






z b a 


/ 

(96) 


Method b 


Ref c 


UVB 


C 


3.1 x 10 4 


2.3 x 10" 


-7 


3 x 10 5 


2.2 x 10- 


-6 


0.0026 


8 


Civ HPOD 


(1) 


soft 


C 


4.9 x 10 4 


3.6 x 10- 


-7 


4.7 x 10 5 


3.4 x 10- 


-6 


0.0041 


13 


C iv HPOD 


(2) 


hard 


Si 


4.6 x 10 4 


3.4 x 10- 


-7 


1.2 x 10 6 


8.6 x 10- 


-6 


0.0100 


30 


Si iv CPOD 


(3) 


soft 





2.4 x 10 5 


2.0 x 10- 


■6 


4.6 x 10 5 


3.4 x 10- 


-6 


0.0040 


12 


Ovi 


(2) 


hard 





7.0 x 10 5 


5.0 x 10- 


■6 


1.2 x 10 6 


8.4 x 10" 


-6 


0.010 


30 


Ovi 


(2) 


soft 





1.7 x 10 5 


2.3 x 10- 


-6 


5.2 x 10 5 


3.9 x 10" 


-6 


0.0046 


13 


Ovi 


(4) 


hard 


Summary 


3.5 x 10 5 


2.5 x 10- 


-6 


5.8 x 10 5 


4.2 x 10- 


-6 


<0.0050 


< 15-30 









a At z = 2, all the baryons are in the forest, and the mean baryonic metallicity Z b corresponds the mean metallicity of the IGM (by 

mass). 

'HPOD: pixel optical depth using Hi; CPOD: pixel optical depth using Civ. 
c fl1 ISchave et al.l (1200311 . (21 ISimcoe et al.1 <2004) . (3l lAguirre et al.1 <2004) . f4) lBergeron fc Herbert-Fortl ll2005h . 



Z b ~ 0.004, 

about 12% of the metal budget. Using a 'soft' UVB gives 
[O/C] ~ 0.5 which is more in line with the relative abun- 
dance s seen in other me tal-poor environments such as halo 
stars l|Cavrel et al.ll2003) . and increases the median [O/H] by 
~ 0.15 dex, but a calculation of Qo u nder this assumption 
was not provided. If the [C/H] values of lSchave et al.1 (|2003l ) 
are used with [O/C]=0.5 (which is consistent with results 
using the pixel optical depth technique; see Dow-Hygelund 
et al., in prep.), we would obtain Z b ~ 0.01, or ~ 30% of the 
metal budget. 

In stead of looking at [O/H], iBergeron fc Herbert- Fort 
l|2005h searched for O VI where the identification of the sys- 
tems is done with Civ, i.e. independently of Afjjj. These 
studies found two populations of O vi absorbers. The first 
population is metal poor (with [0/H]~ —2.0) and has 
narrow line width (b < 12 kms -1 ), indicative of pho- 
toionization. The second (and new) population is much 
more metal rich with [0/H]~ —0.33, and has larger line 
widths. Globally, they found that Qq vi = 3.5 x 10~ 7 , 
of which the metal-rich population contributes 35%. Us- 
ing a ionization correction Ovi/O=0.15 (assuming a hard 



UVB). IBergeron fc Herbert-Fortl |2005l l infer a Oxygen den- 
sity Qo = 2.3 x 10 -6 , which corresponds to: 

pz.forcst ^ 5.2 x 10 5 -^- M Mpc" 3 , (12) 
0.60 

Z b ~ 0.0046, 

or 13 % of t he metal budget, similar to the estimate of 
iPettinil (|2006l ). 

The results from the literature are summarized in Ta- 
ble [2] The Lya forest mean metallicity (by mass) is Z b ~ 
0.005-0.010 (depending on the UVB model assumed and 
tracer element used) and indicates tha t it holds ~ 1 5 — 30% 
of the z — 2 metal budget (see also |Pettinill2006l ). Using 
carbon as a tracer leads to a somewhat smaller estimate . 

Assuming that the intergalactic metal budget is domi- 
nated by warm (~ 10 4 K) photoionized gas, the contribution 
of the forest is 15 (30)% depending on the UVB. If the UVB 
were a bit harder (or softer), it would decrease (increase) 
slightly 1 or 2 elements, but the other element (s) would in- 
crease (decrease). In particular, in order to have a UVB that 
yields a [Si/O] ratio consistent with type II SNe, there is lit- 
tle room to change the 15-30% contribution. 

However, it is important to note that if a significant 
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reservoir of metals is hidden in hot (> 10 5 K) collision- 
ally ionized gas, these would evade detection in C IV and 
Si IV, so that the ionization corrections employed in the cal- 
culations cited above would underestimate the true metal 
content. For example, in sim ulations including feeback by 
lOppenheimer fc Pavel (|2006l ) , heating of the gas hides a sig- 
nificant fraction of the carbo n mass, so that the ionization 
fraction of carbon is than in ISchave et all (|2003l ) by a fac- 
tor of almost three. This might bring the carbon metallicity 
more in line with silicon and oxygen; which might be af- 
fected by prevalent hot gas to a lesser (though currently 
unquantified) degree. 



2.4.2 Lyman Limit Systems (LLS) 

Lyman limit systems (LLS) with Ajjj > 10 17 cm -2 are 
prime candidates for harboring the missing metals. Indeed, 
some of them are highly ionized (because they have a lower 
Hi column density and might therefore not be sufficiently 
self-shielded), and i f a small fraction of LLSs are metal rich 
(as already seen in ICharlton et all 120031 ; iDing et alj 120031 ; 
iMasiero et al.l [20051 ). t hey could contribute significantly to 
the m etal budget (e.g. IPeroux et al.ll2006l ; IProchaska et al.l 
l200d ). 

However, no model-independent constraints exist for 
LLS with 10 17 < iV Hl < 10 19 cm , but progress is being 
made for absorbers with 10 19 < iVjjj < 10 20 ' 3 cm -2 , also 
called sub-DLAs, since dampi ng wings are clearly visible at 
such HI column densities (e.g. IPeroux et al.l 2006). 

Based on direct measures of the neutral gas mass and 
metallicity of a sample of sub-DLAs, iKulkarni et al.1 (|2006l ) 
calculated the amount of metals of these systems to be Z\, = 
6.4xl0~ 4 ,orfi z ~5.1xl0~ 7 using no ionization correction 
/, i.e. / = 1.0 where / = 1/(1 — x) for an ionization fraction 
x. This corresponds to: 



Pz.sub-DLAs > 6.9 x 10 4 ( 



1.07 0.15Z, 







M Mpc" 



(13) 



or about 2% of the metal bud get. Eq. | 13 | is a lower li mit since 
x — was assumed. However . iKulkarni et al.l (|2006T > empha- 
size that x is model dependent and that different sub-DLAs 
might ha ve very different ionization f ractio n even at similar 
iVjj : (see iDessauges-Zavadskv et all (120061 ) for a sample of 

13 ionized fr action estimates) . 

Indeed, IProchaska et alj (|2006l ) deduce x =0.9 for one 
of their sub-DLA for which they run photo-ionisation mod- 
elling and assume £=0.1 for the other one. Based on this 
single measure, they estimated the amount of metals in sub- 
DLAs (i.e. assuming x=0.9) to be f2z ~ 5 10 -6 or: 



2.5 Damped Lya absorbers (DLAs) 

Because (i) S7jj l appears t o evolves by at most a fac- 
tor of two to redshift zero (|Rosenberg fc Schneideil 120031 : 
IZwaan et al.ll2005h. and that (ii) the metallicity evolution 
in DLAs ilProchaska fc Wolfe! I2OO0I: IKulkarni fc Falll 12051 : 
IKulkarni et all 120051 : IProchaska et al.ll2003l ) is much milder 
than the metallicity evolution seen in galaxies (Lilly et al.l 
2003), one is forced to conclude that DLAs 'do not trace ev- 
erything', but are just tracing the gas in the same physical 
conditions at all redshifts. 

Regardless of the origin of the gas, from an obser- 
vational standpoint, DLAs are twenty times more metal 
rich than the forest fPettini et al.lll999l: IProchaska fc Wolfel 
Il997l : IVladilo et alfeod : IProchaska et al.ll2003h. But given 



than they account for 2-3% of the baryons 1 Peroux et al.l 
I2003D with fi D LA = 1.5 x 10 8 /!.7 O M Mpc" 3 , only 5% the 
metals, i.e. 



Pz.dla 



2 x 1O 5 M Mpc -3 



(15) 



/'Z.-u'. ] iLAs < 6.8 X 10 J ( — 



0.1Z G 



M Mpc~ 



(14) 



are i n DLAs given their mean metallicity of — 1.15 or O.O7Z 
(e.g. iPettini et ailll999l : IKulkarni et"alll2005l ). 

There seem to be little amount of dust in DLAs 
as measured directly from depletion pattern (e.g. 
IDessauges-Zavadskv et al.l 20021) or from the attenuation 
of the QSO light (e.g . iMurphv fc Liskel 12004 lYork et al.1 
120061 : 1 Wild et al.ll2006l) . The latter method indicates that 
E(B — B) < 0.01 in DLAs. However, they have been claims 
that dusty DLAs would be missed from optically selected 
quasa r samples altogether as argued by IVladilo fc Perouxl 
(2005). These authors find that the missing fraction of 
DLAs is a strong function of the limiting magnitude of the 
quasar sample and they concluded that (i) we might be 
missing as much as 30-50% of Ojjj for QSO surveys down 
to r = 20.5, and (ii) averaged DLA metallicities could be 
5 to 6 times higher than currently observed. Up to now, 
the magnitude range r = 19.5 to r = 20.5 has been largely 
unexplored, and it is where most of the bias due to dust 
could be more significant. 

Recently, iHerbert-Fort et all (|2006l ) identified from 
SDSS-DR3 a large sample of 435 metal strong QSO absorp- 
tion line selected using EW(Zn), with properti es similar to 
the metal strong DLA of lProchaska et all (|2003l ). From their 
sample, they infer a global fraction of metal strong DLAs of 
5% down to r = 19.5. 

In order to estimate the additional contribution to pz 
from metal-rich dusty DLAs . we examine the v arious pos- 
sibilities. In the scenario of IVladilo fc Perouxl (|2005l ). the 
averaged metallicity of DLAs should be 5 times higher than 
the observed mean DLA metallicity of Z/Zq ~ 0.07. Given 
their estimate of a 'dusty DLA' fraction missed in current 
surveys of 30% (by mass) , the metallicity of this populations 
might be close to be around solar (but not much higher). If 
we take their metallicity to be Z ~ O.5Z , 7 times that of 
'traditional' DLAs, the amount of metals in this population 



or 17% of the metal budget. This result used a mean metal- 
licity for the entire population of 0.1 Zq from IPeroux et all 
(|2003t ). 

Given the potentially large reservoir of metals in 
sub-DLAs, progress in this field is adva ncing rapidly 
jPeroux fc et al.ll2006l : IProchaska et al.ll2006l ). 



6 iFox et al. discussed the contribution of the hot gas in 

DLAs probed by Ovi to the metal budget. This amounts to > 
0.7% of the metal budget and could be higher depending on the 
ionization correction. 
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Table 3. Metals in the IGM. 



Class 


Pz 

(M Q Mpc" 3 ) 


Hz 




z b 

(Zq) 


PZl 'pZ. tot 3 

(%) 


Ref b 


Note 




DLA (JV Hl > 10 20 - 3 ) 


2.0 x 10 5 


1.5 x 10" 


(i 


0.0018 


5.0 


1 






sub-DLAs (10 19 < N Ul < 10 20 3 ) 


6.8 x 10 5 


5.0 x 10- 


-6 


0.0060 


<n 


2 


with x 


= 0.9 


sub-DLAs (10 19 <N Ul < 10 20 3 ) 


6.9 x 10 4 


5.1 x 10" 


7 


0.0006 


>2 


3 


with x 


= 0.0 


LLS (10 17 < N Hl < 10 19 ) 


















Forest (JV Hl < 10 17 ) 


5.8 x 10 5 


4.2 x 10" 


-6 


0.005 


<15 


4 






IGM Metals: Total 


1.5 x 10 4 


1.1 x 10" 


5 


0.011 


< 37 


excluding LLS 







"References: (1) Eg. US] (21 Prochas ka et al.l l|2006h . (31 iKulkarni et all j2006t) . (4) Tabled 



of dusty DLAs is estimated to be: 

PZ, DLA+ = 30% ■ Q.DLA ■ 0.5Zq 

~ 4.7 x 10 5 



0.30 0.5 ^Mpc" 3 ^) 

or about 12% of the metal budget; a non-negligible contri- 
bution. If the fraction of m etal rich systems is closer to ~5% 
IIHerbert-Fort et al . 2006) with a mean metallicity of 0.3 Zq, 
then Eq.[l6]is reduced by a factor of 10, and the contribution 
of dusty DLAs is not significant. This population of dusty 
DLAs could contribute significantly to the metal budget. 



In Tabled we summarize the metal budget for the IGM. 
The sum for DLAs, sub-DLAs and the forest reaches about 
< 35 to <40%, which is barely what is required to close the 
metal budget. 

2.6 QSOs and AGN feedback 

The contribution of QSOs to the cosmic metal budget is 
usually not included due t o their very low nu mber density, 
nqso < 10" 6 Mpc" 3 (e.g. ICroom et alj|2004h. But. QSOs 
appear to h ave reached solar metallicity (e.g. lllamann et"al] 
2002, l2004h based on independent analyses of quasar broad 
emission lines and intr insic narrow absorption lines (see also 
iD'Odorico et al .120041 1. Given that their life time is relatively 
short (r < 10 Myr, i.e. duty cycle large > 100), their con- 
tributions can be much larger if they contain large amount 
of gas. From CO observations, it has been note d that QSOs 
conta in less gas than SMGs by a factor of ~ 4 (|Greve et al.l 
2005), so we expect their contribution to be less than that 
of the SMGs. Typic ally, QSOs have ~ 10 10 M of gas (e.g. 
Hainl ine et al.ll2004h . We find that, the comoving amount of 
metals in QSOs is then: 

Pz,qso = 2 x lO^X^i Mq Mpc" 3 (17) 

or less than 1% of the expected metals. 

Thus, even if QSOs expel most of their gas, they can- 
not contribut e to the metal budget a t more than the percent 
level. In fact. iNesvadba et alj l|2006l l studied the kinematics 
of the outflow powered by the radio galaxy MRC1138— 262, 
and estimated an outflow rate 300-400 Mq yr _1 . From the 
outflow rate, they estimated the contribution to pz from 
AGN feedback assuming that all QSOs undergo a power- 
full radio phase similar to that of MRC1138— 262 (using a 



duty cycle of 300). They found that they could contribute 
to 0.1— 30xl0 4 M Mpc" 3 of metals, i.e. at most 10%, and 
the most likely value is 1 x 10 4 M Mpc -3 or 0.3% given 
that QSOs outflows are not likely to expel more metals than 
their gas reservoir (Eq. I17[) . We thus view that AGN feed- 
back is unlikely to contribute significantly to the enrichment 
of the IGM, but wi ll significantly affect t he evolution of 
mass i ve objects (e.g. | Croton fc et al. JI2005I : INesvadba et al.l 
120061 : iBest et al.ll2006l l 



3 HOW MANY METALS ARE EJECTED 
FROM SMALL GALAXIES? 

In this section, we use simple analytical calculations based 
on the effective yield y c g = Z/\n(l/p) (Eq. [Q) to argue that, 
at z = 2., 25-50% of the metals are 'outside'[j galaxies. First, 
we test our method at z = in section 13.21 and then repeat 
it at z = 2.5 in section EP1 



3.1 The contribution of metals lost from galaxies 

At 2 = 0, evidence for enriched material being ejected 
by galaxies comes fr om the mass (or luminosity)- 
metallicity relation ( e.g. Garnctt 2002; Pil vugin et al. 1 12004 
ITremonti e t al. 2004) . Essentially any chemical model would 
predict that a galaxy (as it turns its g as into star) reaches 
solar metallicity (e.g. Edmund's! Il990l l even with infall of 
primordial gas. Once the energy of SN is larger than the 
gravitational energy of the gas, the remaining gas could be 
expelled, resulting in a mass-metallicity rela tion since th e 
gravitational energy will depend on the mass l|Larsonlll974l ). 

Another signature of gas losses come from the effective 
yield (Eq. [TJ, which measures how far a galaxy is from a 
'closed-box' evolution. A galaxy that has evolved as a closed- 
box would obey a simple linear relationship between gas 
metallicity and gas fraction, the slope being the effective 
yield, i.e. Z = y c g x ln(l//u) (Eq.[T}. 

Fig. [2] shows the effective yield as a function of rota- 
tional vel ocity V c . The so lid line shows the relation parame- 
terized bv lGarnettlll2002ll. and the dashed line shows the re- 
lation obtained by ITremonti et all (|2004T l using SDSS data. 
The shaded area represent the allowed range from the ob- 
servations. 



7 By this we mean not in stars, and not in the ISM. The metals 
could be in the halos of galaxies or in the IGM proper. 
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Figure 2. Effective yield (y c g) normalized to the value expected 
from a close b ox model yp. T he solid line shows the relation pa- 
rameterized bv lGarnetd l2002h, and the dashed line shows the re- 
lation obtained bv lTremonti et al] J2004T) using SDSS data. The 
shaded area represent the allowed range from the observations. 



This effective yield has two important properties rel- 
evant for t his paper . Firs t , as it has been shown many 
times (e.g. lEdmundsl Il990l ; iDalcantonlkoOrJ ). outflows can 
lower j/eff effectively. Second, how much the effective yield 
has departed from the closed-box expectation (given b; 
true yield) is (e.g. iGarnettl |2002|; iTremonti et al.l 1200 



iPilvugin et~ai1 12004 iBrooks et all l200fj ) a measure of the 
minimum amount of metals that were lost in the last star- 
burs t driven outflow, since, tech nically, inflows can increase 
2/eff (|Koppen Sz Edmunds! Il999h after an outflow episode. 

It can be argued that the ratio between the effective 
yield y B a and the true yield yo is a measure of the lost metals 
by approximating the effective yield when the gas fraction 
H = M gas /(M gas + Mi,) is close to 1, in which case: 



Vett 



M z 



(18) 



where Mz is the mass of metals in the gas (= Z M gas ). One 
sees that the effective yield is independent of how much gas 
is in the system, but any process that removes metals (as 
in metal-rich outflows) will reduce y B s in proportions to the 
metals lost since Eq. [TS]is linearly proportional to the mass 
of metals in the gas phase. 

One can show more directly that the ratio y e s/yo is 
a measure of the mass fraction of the lost metals (see Ap- 
pendix . A galaxy that experienced an outflow episode 
has a remaining mass of metals M%- The mass of metals 
lost (Afi? st ) with respect to that of the closed box evolution 



(Mz ), normalized to the mass of metals left M$ is (Eq.fM 



Mi 



M% 



Ml 



Ml 



yo AC 

^ff Af» r 



(19) 



where M£ ar and M^ r are the baryonic masses for the closed 
box evolution and for the wind scenario, respectively. 

In general, the baryonic mass is not affect by outflows, 



1, and 



Ml 



yo 



(20) 



From Eq. 1201 one sees that the ratio between yo and y c s 
gives the fraction of metals that were lost from the system. 
For instance, if j/off/yo =1/2 then the amount of metals lost 
is equal to those remaining, i.e. 50% of the metals produced 
were lost. 

This property of j/ e fi can be used to compute the mass 
of metals eject ed from galaxies of a given luminosity Lb (e.g. 
IGarnettl liool ) as follows: 



Mo.iost = 12- (L B ) Lb^- 
ti Lb 



yo 



M . 



(21) 



where 0/H(Lb) is the luminosity-metallicity (L-Z) relation, 
M^ I Lb the stellar mass ratio, the factor 12 converts the 
oxygen ratio to a mass ratio taking into account the He 
fraction, and y e g the effective yield. 

To illustrate Eq. 1211 we compare its prediction to the 
obser vations of NGC1569 (|Martin. Kobulnickv. Heckmanl 
2002), where the amount of Oxygen in the outflow was mea- 
sured. This nea rby dwarf ga laxy has a baryonic mass of 
~ 2 x 10 8 M Q ijMartirj 11998 ). a metallicity of 0.2Z Q , and 
entered a starburst phase 10-20 Myr ago. The cumulative 
effect of ~30,000 supernovae in the central region is driving 
a bipolar outfl ow seen exten ding to either side of the disk in 
Ha emission ([Martin 19981) and in X -rays. Using Chandra, 
iMartin. Kobulnickv. Heckmanl (|2002r ) measured the amount 
of metals in the outflow from the X-ray emitting gas, and 
found that the hot wind carries 3.5 x 10 6 Mq of gas which 
includes 3-5 x 10 4 Mq of Oxygen. The gas and the cur- 
rent stellar masses are both uncertain in this galaxy. For 
a gas fraction of ~40, 50, & 60% and given its baryonic 
mass, the effective yield is about 0.4, 0.5, and 0.7, respec- 
tively. The corresponding amount of Oxygen in the ISM is 
1.5-3 x 10 5 M Q . Using Eq. [20] the mass of Oxygen lost is 
about 4, 2.5, and 1. 2 xlO 5 M R . This is larger than the d i- 
rect measurement of IMartin. Kobulnickv. Heckmanl (|2002f ). 
However, the yield reflects the cumulativ e effect of all the 
past bursts. In fact. lAngeretti et al.l (|2005l ) showed that this 
galaxy likely experienced three separate starburst phases, 
giving a mass loss per starburst phase of 1.3 to 0.4 x 10 5 Mq. 
This is not far from the amount of metals in th e curre nt out- 
flow found by M artin. Kobulnickv. Heckma n (2002) given 
that some additional amount of oxygen can be in a much 
hotter (T >> 10 6 K) gas difficult to detect. 

Eq. [21] can be used in combination with a luminosity 
function (LF) to compute the global amount of oxygen (per 
unit volume per unit magnitude) ejected into the IGM as a 
function of luminosity : 

dMo 



dVdmag 



= M OAost (M B )4>(M E 



(22) 



whose integral with respect to magnitude gives the total co- 
moving density of oxygen lost po from the ISM of galaxies. 

In section I3~2l we perform the integral (whoe input de- 
tails are presented in the AppendixfX]) at z — using either 
the Lb luminosity function or the stellar mass Af* function. 
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Wc find that both approaches give very similar results. Fur- 
thermore, we will show that the amount of metals lost from 
galaxies corresponds to what is observed in the intra-group 
and cluster medium and conclude that our methodology is 
sound. 



3.2 Metal loss at z = 

The solid line in Fig. [3] (left) shows the amount of oxygen 
lost (i.e. ejected in their halos, and/or in the IGM) in units 
of Mq Mpc -3 produced by z — galaxies computed using 
Eq. [22] This solid curve peaks at ~ 0.2L* and shows that 
the bulk of the metals (given by the median or the peak) is 
ejected by sub-L* galaxies. The dotted line show the same 
but per unit lumino sity (right axis ) instead of magnitude, 
for comparison with iGarnettl (|2002l ). 

The solid line in Fig. [3] (right) shows the cumulative 
distribution of po.lost using Garnettl (|2002l Vs parameteriza- 
tion (Eq. IB1|I of j/eff- We also show the 50th (90th) per- 
centile as open triangle (pentagon). The dashed curve shows 
the result of t h e int egration using y e g parameterized by 
iTremonti etafl l|2004h (Eq. [B5J. From this plot, one sees 
that the total amount of oxygen ejected from galaxies is 



Po.ejcctod^ = 0) « 3.0 x 10 6 M Mpc" 3 . (23) 

The oxygen yield is about 60% of all metals l|Samlandl 

1998), so we find that the total co-moving density of metals 
lost from the ISM of galaxies: 



Rejected (z = 0) « 5.0 x 1Q 6 M Mpc~ 



(24) 



This correspond to ~25% of the z = metal budget (Eq(4j) 
or to Zt = 0.045. The ingredients that went into Eq. 1241 are 
presented in Appendix IB II 

One can perform the same calculations directly in terms 
of stellar masses M» using stella r mass functions. Using the 
results of ITremonti et all l|2004h for y eS (Eq. IB5|I and for 
the mass-metallicity (M*-Z) relation (Eq.lB4|). the baryonic 
Tully-Fischer relation fro m dBell fc de Jondl200lh and the 
stellar mass function from lRead fc Trenthaml (|2005r ). we find 
that the comoving oxygen density is po, ejected (z = 0) ~ 
2.0 x 10 6 Mq Mpc -3 , corresponding to a metal density of 



PZ,ejected(2 = 0) 



3.3 x 10 b M Mpc' 



(25) 



This is close to our estimate based on the LF (Eq. I24|l and 
very different assumptions. The ingredients that went into 
this estimate are presented in Appendix IB2I 

Direct measurements of the plasma associated with 
galaxies (outside rich clusters) are still quite uncertain, but 
their contribution to the baryon budget ap pear to be sig- 
nificant. According to Fukug ita et al.l (|l998D. the mass den- 
sity of the group plasma is low f2 gr = 0.005 for h — 0.70. 
The metallicity of this plasma can only be measured in the 
brightest X-ray groups and ra nges from 0.1 sola r to 0.6 so- 
lar. According to the review of lMulchaevI |2000), this intra- 
group medium has a mean m etallicity similar to that of clus- 
ters, i.e. 1/3 solar (see also iFinoguenov et al.l [20061 ) . As a 
consequence, the amount of metals observed in groups is, 
pz, g rou ps (z = 0) ~ 4.3 x 10 6 M Mpc -3 , around 18% of the 
z — metals produced by type II supernovae. 

The contribution of clusters is about half that of groups 
since the cluster mass density O c i is about half that of groups 



l|Fukugita et al.|[l998T l and both groups and cluster plasmas 
have similar metallicities. 

All in all, the amount of metals outside galaxies ob- 
served in the intra-group and intra-cluster medium is 



PZ,plasma(z = 0) ~ 6 X 1Q 6 M M P C" 



(26) 



around 25-30% of the z = metals, close to our estimate 
based on the effective yield (Eq. [24}. This leads us to con- 
clude that our procedure to estimate the amount of metals 
outside galaxies (i.e. not in stars and not in the ISM) is sen- 
sible, given the uncertainties. The agreement is in fact re- 
markable, since the normalization of Eq.[22]d e P e nds strongly 
on the normalization of the LF, and on the effective yield 
method. 



3.3 Metal lost at z = 2.5 

We now turn towards the redshift of interest z ~ 2.5, and 
perform a similar calculation as in section 13.21 The details 
of our assumptions are ex plained in secti on IB3I Briefly, 
using the j/ e ff relat ion from IGarnettl (|2002f l (Eq. [BlJ, the 
L—Z relation from lErb et all (|2006h (0.3 dex offset com- 
pared to z = 0), the -B-b and luminosity function from 
ISawicki fc Thompson! ([20061 ) ■ the B-band Tully-Fischer rela- 
tion (TFR) offset by 1 mag (see section [B3j 0, and a M/L B 
ratio equals to 2, we find that 



PO,ejectcd(z — 2.2) 
PZ,ejectcd(z = 2.2) 



1.3 x 10 6 M Mpc -3 , (27) 
2.1 x 10 6 M Mpc" 3 , (28) 

which corresponds about ~ 50% of the expected metals 
(Eq. O and to a mean metallicity (if we spread these metals 
over all the baryons): 



Z b = 0.018. 



(29) 



The solid line in Fig. [4] (left) shows the amount of oxy- 
gen lost in units of M Mpc -3 produced at z = 2.5 com- 
puted similarly as in Fig. [3] The solid curve peaks at ~ L* 
and shows that the bulk of the metals (given by the median 
or the peak) is ejected by sub-L* galaxies. The dotted line 
show the same but per unit lumin osity (right axi s) instead 
of magnitude, for comparison with IGarnettl (|2002T l . 

The solid line in Fig. fright) shows the cumula- 
tive distribution of po,ioat- The 50th (90th) percentile are 
shown as open triangle (pentagon) . The dashed curve shows 
the result of the int egration using y e g parameterized by 
ITremonti et~aU (|2004l ) fEq. lBB) . 

These plots indicate that sub-L* galaxies have ejected 
enough metals (50%) to close the metal budget, some (1/3) 
of which is already detected in the forest according to our 
results in section [2.4.11 and the remainder likely being in a 
hotter phase. In particular, low mass (L B < ±L*{z = 2)) 
galaxies are responsible for 90% the production of these hot 
metals, the median being at L ~ jqL* b (z = 2). Given that, 
at z = , the mass outflow rate is usually comparable to 
the SFR llLehnert fc Heckmanlll996bl ; Irleckman etalll2000l : 
IVeilleux et al. 2005, for a recent review), an d that direct 
evidence for winds for LBGs are numerous fe.g lPettini et all 



8 no-evolution of the B-band TFR decreases our results by a 
factor of ps 2. 
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Figure 3. Left: the solid line (and the left y scale) show the amount of metals lost as a function of Mb using our 2 = effective yield 
calculation and Eq. 1221 The dashed line (and the right y scale) show the relative amount of metals lost per unit luminosity. The solid 
line shows that the metals are most likely coming out ~ L* B galaxies. Right: The solid curve shows the cumulative integral of the solid 
curve shown in the left panel. The dashed curve show the result if one uses ^ an( j y cS parameterized by iTremonti et al ] <|2004h . The 
50th and 90th percentiles are shown with an open triangle and pentagon, respectively. The 50th percentile is at ~ O.IL^ galaxies and 
shows that dwarf galaxies eject more metals into the IGM than > L* galaxies. Percentages of the z = metal budget (Eq. [4J are shown 
by the horizontal dotted lines. This plot shows that the metals ejected from galaxies account for 20-30% of the metal budget from our 
effective yield calculations. Since this amount is very close to current estimates of the observed metal mass in plasmas (see text), we 
conclude that the y e g methodology is sound. 



log [L„/L-] 
-2 
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log [L/L-] 

Figure 4. Same as Fig.|3]but at z ~ 2.5. Left: the solid line (and the left y scale) show the relative amount of metals lost (Mq Mpc — 3 ) 
as a function of luminosity using our z = 2.5 effective yield calculation (see text) and Eq. 1221 Right: The solid curve shows the 
cumulative integral of the dotted curve shown in the left panel. The dashed curve show the result if one uses |j and y c g parameterized by 
ITremonti et al. | d2004h . The 50th and 90th percentiles are shown with an open triangle and pentagon, respectively. The 50th percentile 
is at ~ 0.1L|j galaxies and shows that dwarf galaxies eject more metals into the IGM than > L* galaxies. Percentages of the z = 2.5 
metal budget (Eq. [5j are shown by the horizontal dotted lines. This plot shows that about 50% of the metals have been ejected from 
galaxies, 90% of which by < L* galaxies, which is enough to close the metal budget: > 40% are missing when one accounts for z = 2.5 
galaxies (Table [TJ. 



l200ll ; lAdelberger et all2003l : IShaplev et al.ll20oj) . our result 
may not be very surprising. 



at earlier times than at z = 0. This limit > 15% is consistent 
with our estimate of ~ 50%. 



Our conclusion that ^50% of the metals have been 
expelled from sma ll galaxies is very consistent with 
ISimcoe et al l (I2OO4T) who treated the universe as the 'ul- 
timate closed box model', and found that the fraction of the 
metals formed in a starburst that are ejected from the galaxy 
into the IGM, / c j is > 15%. They argue that the value of / e j 
is probably higher since gas to stellar ratio are much higher 



4 INSIGHTS FROM SIMULATIONS 

Several approaches have been taken to predict the cos- 
mic metal abundances and its evolution analytically. Early 
chemi cal models trace only the cosmic mean metallicity Z 
(e.g. iPei fc Fall|[l995l ; lEdmunds fc Phillipps|[l997l : |Pei et al.l 
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1999). Rece nt chemical models (e.g. ICalura fc Matteuccil 
2004, 2006) have been able to keep track of each ele- 
ments over cosmic times, but in typically only 3 types 
of galaxies, spheroids, spir als and dwarfs. In the study of 
ICalura fc Matteuccil l|2006l ). where they modeled spheroids 
and dwarfs only, spheroid formation at z ~ 5 produces most 
the IGM metals, such that the IGM metallicity is fairly flat 
from z = 4.5 to the present. 

Many numerical models have been developed in or- 
der to simulate the distribution of the metals in the 
IG M. But, even in models that do include galactic winds 
(e.g. l Cen fc Ostrikerll999l:lAguirre et al.l200ll ; lTheuns et all 
12003 : iBertone et alj 120051 ; ICen et al.ll2005l ). the predicted 
IGM metallicity is too high compared to the mean metallic- 
ity of the forest, whi ch is very low (0.005 Zq). For instance, 
IBertone et ail (|2005l ) modeled galactic winds analytically in 
N-body simulations and found that winds should have en- 



riched the IGM to a metallicity of 10-10" 



about 10 



times larg e r than what is observed, a conclusion reached by 
Cen et all (|2005l) using very different simulations. Recently, 



Bertone et al. I (|2007l ) investigated the fraction of metals lost 



from galaxies in galaxy formation models. They specifically 
looked the amount of metals that are permanently lost from 
the host galaxy, and computed the distribution of metals lost 
as a function of virial mass M v i r . They found that at z = 0, 
the 50th precentile of this distribution is at ~ 10 11 ' 5 Mq, 
or about 0.8dex smaller than L* (using M* ir — 10 12 ' 3 Mq). 
This is rather close to our estimate shown in Fig. [3] (right) 
of log[L/L*] = 0.6. A preliminary comparison between their 
z — 2 result expressed as a function of luminosity (Bertone, 
S. private communication) and Fig. |4]yields very encourag- 
ing qualitative agreement. 

Keeping track of individual elements in cosmological 
simulations is more difficult. Recently, chemical models of 
individual elements have b e en attached to S PH simulations 
|Samland fc Gerhard! 120031: iKobavashil 120041) ■ bu t are often 
limited to single halos. Kobavashi et alf ~(2006) combined 
chemical evolution models in a full cosmological tree-SPH 
simulations. They find that 20% of all baryons are ejected 
at least once from galaxies into the IGM. Galactic winds 
are found particularly efficient in low mass galaxies. They 
argue that the origin of the mass-metallicity relation is from 
galactic winds. 

Using hydrodynamical simulati ons that incorporate 
metal -enriched kinetic feedback, iDave fc Oppenheimerl 
(2006) concluded that 50% of the z = 2 metals are not in 
galaxies. These simulations use scaling relations that arise in 
mom entum-driven winds (|Murrav. Quataert fc Thompson! 
120051 ) and ionization calculations to calculate the ionization 
fractions of atomic species given the gas density, tempera- 
ture and the ionization field. In order to reproduce the abun- 
danc es of C iv from z — 6 to z — 2 l|Oppenheimer fc Pavel 
2006) generally requires high mass loading factors but low 
velocity winds from early galaxies, in order to eject a sub- 
stantial metal mass into the IGM without overheating it. In 
other words, momentum-driven outflows provide the best 
agreement with observation s of C IV. The simulations of 
IDave fc Oppenheimerl (|2006l ) are in good agreement with 
several of our results. For instance, they find that 20-30% 
of the metals are locked in stars at z = 2, 15-20% are in the 
shocked IGM (akin to Ovi), and 15% reside is cold star- 
forming gas in galaxies (akin to DLAs and sub-DLAs) and 



30% are in the diffuse IGM (akin to the forest). Essentially, 
they find that 40% of the metals are in galaxies, in agree- 
ment with our constraints of > 30 and < 30%. Half of the 
remaining 60% is predicted to be in the diffuse IGM (al- 
though with substantial fraction in hotter gas), and the rest 
are divided between the shocked IGM and the hot halo gas. 



Madau et al. 


(120011): Dave et all ( 


200 ll): IScannapieco et al. 


(2002): Furlanetto & Loebl ( 


2003): iDe Lucia et al. 


(2004): IScannaoieco! (20051): iPorciani & Madaul (20051): 



5 SUMMARY & CONCLUSIONS 

In terms of the z ~ 2.5 metal budget, the dominant contrib- 
utors are the galaxies, the forest and sub-DLAs: 

(i) }t 30% and < 60% of the metals are currently observed 
in z ~ 2.5 galaxies (Eq. [7}, according to our calculations 
in paper I, and paper II extended here and summarized in 
Table □ 

(ii) 15-30% of the metals are in the low column density 
IGM with logNjjj < 17, and < 17% are in sub-DLAs with 
19 < log[7Vjjj] < 20.3 as estimated, in section [5] and sum- 
marized in Table [2] 

(iii) from our effective yield calculations (Fig. [4}, ~50% of 
the metals are predicted to be outside galaxies. This method 
was successfully tested at z = 0, where the census of metals 
outside galaxies is more robust. Furthermore, we find that 
low mass (Lb < \L*(z — 2)) galaxies are responsible for 
90% of these metals. 

(iv) this last result agrees well with observations since we 
can account for ~ 30-45% of the metals outside the ISM 
of galaxies by adding the contribution of the forest and of 
sub-DLAs. We note, that the contribution of absorbers with 
17 < logfAfjjJ < 19 is currently lacking, and that this class 
(plus the sub-DLAs) are potentially a good place to look for 
the remaining metals as current constraints are very limited; 

(v) comparing the amount of metals outside galaxies at 
z = (Eq. [24} and at z = 2.5 (Eq. [28} , one sees that it has 
evolved by a factor of 2. 

Taken at face- value, one could conclude that the metal 
budget is almost closed since 30-60% of the metals are in 
galaxies and < 30-45% are in the IGM. However, given the 
potential overlap between the various galaxy populations, 
we think it is very unlikely that all the metals have been ac- 
counted for. Conservatively, only about >65% of the metals 
have been detected directly, equally spread between galaxies 
and the IGM. 

Our effective yield calculations indicate that metal rich 
outflows from galaxies are likely to be the reservoir of the 
remaining missing metals. A significant fraction (2/3) of this 
is already seen in absorption lines studies and the remaining 
metals are very likely in a hot phase (hotter t han photoion- 
ization tem peratures), a conclusion reach ed bv lFerrara et alj 
(|2005h and bave fc Oppenheimerl l|2006l ). 



6 DISCUSSIONS 

We end this paper by discussing three related questions. 
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• Are the metals seen in the IGM at z = 2 produced by 
z = 2 galaxies? In other words, is the comoving number den- 
sity of 'BX' galaxies large enough and is their star formation 
phase lasting long enough for producing 50% of the z = 2.5 
cosmic metal density? Given the co -moving density of z = 2 
galaxies is n ~ 2 x 10 
their typical SFR are 



h? Mpc 3 jAdelberger et al 



35 M fl yr 1 dShaplev et al 



2004), 



2003), 



ages of 2 x 10 8 yr (jShaplev et al.l 120051 ). and metallicities 
of 0.5 Zq (|Erb et alj|2006h . the co- moving metal density in 
winds from star-forming galaxies is: 



Pz = 



5.3 x 10' 
Z 



G) 



SFR 



0.5Z, 



35M yr" 1 
n 



2 10 8 yr 



2 10- 3 Mpc^ 



M Q Mpc" 



(30) 



where we assumed a w ind mass outflow rate r of 4 times the 
SFR (|Erb et all [20061 ). Eg . 1301 amounts to about ~10% of 
the metal budget, i.e. much smaller than the > 40% missing, 
or the 50% expected to be outside galaxies 0. This is telling 
us that 3-5 bursts of star formation are required per galaxy. 

• A possibility for the remaining missing metals is that 
they are in a phase hotter than photoionization tempera- 
tures, at T > 5 x 10 J K. A natural question is then: 
are there enough energy in the outflows to keep some of 
the gas invisible? The energy for one hydrogen atom (or 
baryon) in the outflow of velocity V ou tfl w is \mH Voutflaw 
The characterisitic o u tflow speed is 400-800 kms -1 (e.g. 
iHeckman et al. 2000; iMartinl I2005D at z = and ~ 
600 kms" 1 at z — 2.5 from the s hift of the the low-ion ization 
lines to the Lya emission lines (|Shaplev et al.ll2003l ). Thus, 



the outflow carries 



( ^outflow y 

Veoo kms -1 / 



keV/baryon of ki- 
netic energy. Whereas, the energy necessary to heat the gas 
to ~ 10 6 K is § kT, which is ~ 120 (j^k) eV per baryon. 
Thus, only 6% (ranging from 3 to 15% for V ou tflow =400- 
800km s -1 ) of the kinetic energy available is necessary to 
heat the gas to T ~ 10 6 K. 

• Are the metals seen in the IGM at z — 2 going to fall 
back or leave the virial radius for ever? If galaxies evolve with 
only outflows at an outflow rate always equal to the star for- 
mation rate, then there would always be as many metals in 
stars as in the IGM at all times, including at z — 0. Thus, 
the ratio between the metal density of the IGM to the metal 
density in stars would have to stay constant. But, at z = 0, 
about 80-90% of all the metals produced by type II and 
type la supernovae are locked in stars, i.e. < 10% of the total 
metal density are in plasmas outside galaxies, i.e. that ratio 
is 1:9. From this work, the metals outside galaxies amount to 
Pz, c jcctcd(z = 2) ~ 2 x 10 6 M Mpc -3 and the stellar metal 
density is /9 Z , st ar S (z = 2) ~ 1-2 x 10 6 M Mpc -3 , i.e. that 
ratio is 2:1 or 1:1. Since this ratio has evolved from 2:1 at 
z = 2to 1:9 at z — 0, we view this as evidence that a signifi- 
cant fraction of the IGM metals have cooled and fallen back 
into gal axies. This conclus i on wa s reached using 2 = SDSS 
data by iKauffmann et all (|2006t ) on very different ground, 
and is consistent with the models of iDave fc Oppenheimen 
(2006). This dramatic change in proportions may be related 



9 We not e that , technically, the outflow 'rate' r inferred by 
lErb et al] ||2006|) gives a true outflow rate of (1 — R) r times 
the SFR, where R is the mass fraction in massive stars (~ 0.3). 



to a shift in the dominant phase of star formation: where, 
at z = 2.5, bursty (i.e. producing outflows) star formation 
dominates due to rapid accretion and/or merger, while more 
quiescent mode of star formation dominates at smaller red- 
shifts due to more slower accretion. This is very reminiscent 
of the redshift dependence of the two modes of gas accre- 
tion discussed in ?. At high redshifts, cold accretion along 
filaments dominates, while at low-2, the hot accreation dom- 
inates. Both modes qualitatively reproduce the rapid vs. qui- 
escent mode of star formation. 
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APPENDIX A: EFFECTIVE YIELD 
CALCULATIONS 

In Appendix I All we justify Eq. 1201 which is central to this 
paper and describe the behavior of Eqs. [21] and [22] in Ap- 
pendix [Al 

Al Fraction of metals lost 

One can show that the fraction of metals lost (with respect 
to a closed box model) is proportional to yo/y c g — l. Consider 
a galaxy that has evolved for a given amount of time, untill 
some gas fraction /i remains, and experienced an super-nova 
driven outflow. If it had evolved as closed box instead (untill 
the same gas fraction), the amount of metals it would have 
is given by 

Mf = Z ch (Mg M +M°L a ), 

= j/oln(l/M)(Mis b M + M s t rs ), (Al) 

where Mf SM + M^ als is the baryonic mass M^ T and j/o is 
the true yield. 

Instead, it is observed to be metal poor for its gas frac- 
tion /x. The amount of stars M sta rs with or without outflow 
would be the same as it depends only on the star formation 
rate. In the case of outflows made purely of SN ejecta, the 
amount of gas Mism left after the outflow episode is also 
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the same. This SN ejecta simply left the ISM and failed to 
enrich the remaining gas and stars. If the outflow entrained 
some fraction of the ISM, the gas mass may differ from that 
of the closed box evolution. 

In the most general case, the total amount of metals 
after an outflow episode is: 



= y a s ln(l/ n) (M^ M + M s ™ 

ars / 5 

where y e g is the effective yield. 

The amount of metals that were lost is M 



(A2) 



lost 
Z 



= M% 



. The fraction of metals lost relative to what is left MS 



is 



MS 



Ml 



My 



Ml 
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bar 
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(A4) 



Since the baryonic mass is not changed by outflows made of 
100% ejecta, and is likely to be close to its original baryonic 
mass even in the case of outflows made of entrained material, 
Eq. IA3I gives the ratio between the amount of metals lost 
from galaxies to that the amount of metals left. 



A2 Dependence on the parameters 

The results presented in Figures [3]-??, can be understood 
from Eqs. (2D and 1221 The core ingredient is the effective 
yield, y a g. Because y e g increases with mass or luminosity 

the term — 1 decreases with Lb- In the case of Eq. IB 1 1 

(as it can be estimated from Fig. [5]) the effective yield 
goes as y c g oc L^ 1 ' 3 . In the case of Tremonti's parame- 
terization (Eq. IB5[) . using Vrot oc L 
from the B-band Tully-Fischer relation. In what follows, 
we parameterize y c g oc L%,. Because the LZ relation in- 



creases with Lb (e.g iLequeux et al.l 11 979; Skil lman et al.l 
Il989l ; IZaritskv et al.lll994l ; iPilvugin et al ll2004h . the factor 
§L B in Eq.l2T1is a steep function of Lb, it is oc L B +2 5x0 16 oc 
L B A - As a consequence, the amount of metals lost per galaxy 

" -1.4. 



according to Eq. [21] goes as Mo, lost oc L 13 

Afo,iost(L B ) oc L+°- 73 for /3 = 2/3, 
Afo.iostOM oc Lj 106 for p = 1/3, 



r t 

B B 



(A5) 
(A6) 



which means that a bigger galaxy ejects more metals than 
a dwarf in absolute terms. 

The shape of the comoving density of metals per unit 
luminosity (given by Eq. I22[) can now be estimated since 
it is given by the amount of metals per galaxy times the 
luminosity function. At the faint end of the LF, a ~ —1.0, 
such that the comoving density of metals per unit luminosity 
will be almost constant: 



oc L~°' 2r for = 2/3, 



dMp 
dVdL 

^L « ^°W = V3. 



(A7) 
(A8) 



Table Bl. Parameters used in section [3.21 and Fig. [3] 



Parameter 



Value Error Relative Error 
Error to pQ 
(1-a) (%) (l-<x) 









9 c +1.3E6 

Zb -l.iEe 


LF 0* 


0.0088 


0.002 


+9.1E5 
-9.1E5 


LF 


-20.27 


0.1 


+1.7E5 
-1.7E5 


LF a 


-1.05 


0.05 


+3.0E5 
-2.7E5 


TF slope 


-7.3 


0.5 


+3.4E5 
-2.8E5 


TF intercept 


-20.11 


0.15 


+5.0E5 
-4.6E5 


LZ slope 


-0.16 


0.01 


+2.0E6 
-1.3E6 


LZ intercept 


-6.4 


0.10 


+9.9E5 
-7.9E5 


M/L 


6 


2 


+1.3E6 
-1.3E6 


Uncertainty to Eq. 1241 






+3.1E6 
-2.5E6 



increase steeply: 

dMo r ~o.7 f o 0/0 

oc L B for p — 2/3, 



d\/dmag 
dM 



oc Lb for P = 2/3. 



(A9) 
(A10) 



dFdmag 

The steep exponential cut off of the bright end of the LF 
will make the contribution of large galaxies to the comoving 
density of metals per unit magnitude drop sharply. Thus, 
somewhere between L << L* and L >> L*, the distribution 
dvdmag C^l- EH wm peak at an absolute magnitude Mb < 
L*. This conclusions does not change even ifM/Ls depends 
slightly on Lb- 



APPENDIX B: ASSUMPTIONS AND 
INGREDIENTS 

In this Appendix, we detail the ingredients used in our calcu- 
lations of the global amount of metals ejected from galaxies 
presented in sections 3.2, and 3.3. 

Bl Metals lost at z = 

Here, we list the ingredient that went in Eq. 1211 F irst, for the 
y c ff as a function of rotational velocity Vrot , Gar nettJ (|2002| '1 
parameterized it as (for galaxies with V I0 t < 250 kms -1 ) 

(320 - Kot) 4 



logy C ff = log y - 



0.1 x 10 9 



(Bl) 



0.01122 is the true yield for Oxygen 



where yo = 10 
(O). 

The second central ingredi ent in Eq. 1211 is the well 
known L-Z relation. At z = 0, iGarnettl (|2002l ) character- 
ized it as: 

12 + log -= -0.16 Mb - 6.4 . (B2) 
H 



and the comoving density of metals per unit magnitude will 



Similarly, iTremonti et~aH l|2004l ) found 12 + log § = 
-0.18 Ms - 5.238. 

The third important ingredient is the LF and the 
_B-band TF-relation. We used a LF with a = —1.05, 
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<j> t = 0.017 * If a nd M* - 5 log ft = -19.6 5 
ICroton fc et al.J l|2005h similar to tha t of Liske et al 



from 



2003), 



and the B-band Tully-Fischer (e.g. iTullv fc Pierce] |2000j) 
M B = -7.3 log[2V r ot - 2.5] - 20.11. 

The final and critical ingredient is the value of MILp, F"l 
as it sets the normalization in computing pz, ejected (Eg. [22)1 ■ 
Many ha ve estimated the mean stellar mass-to-light ra - 
tios (e.g. iBell fc de Jond I2OO1I ; iKauffmann fc et al.J 12003). 



The M/L ratio varies strongly as a functi on of B-V color 
Bell fc de Joir3l200ll:lRudnick et al.J l2006). or galaxy type 



Read fc Trenthaml l2005 ) , from 1 to 9. One way to estimate 
(M/Lb) is from the ratio of the stellar density p* to B- 
band luminosity density js ■ There has been many estimates 
of the stellar mass functi o n at z = using various methods. 
For instance, ICole et all l|200il ) computed the global stel- 
lar density p* from the Tf-band luminosity funct ion of 2dF 
;alaxj es, and found p* = 5.6 x 10 s M Mpc~ 3 . lr3ell et all 
2003) used a similar approach (using a 'diet' Salpeter IMF) 
and found that p* = 3.7 x 10 s hro M Mpc" 3 , correspond- 
ing p, = 5.9 x 10 s hro M a Mpc" 3 for a S alpeter IMF, 
very sim ilar to the results ofTCole et all (|200lh . Recent esti- 
mates bv lRead fc Trenthaml (|2005l ) showed that = 0.0028 
(p* = 3.78 10 s M Mpc" 3 ) using a Kroupa IMF, which cor- 
responds to p * = 7 .5 10 s M Mpc" 3 for a Salpeter IMF. 
IPanter et all l|2004[ ) found for a Salpeter IMF ft* = 0.0034 
(p* = 4.6 x 10 s M Mpc" 3 ). We will use 

p* = 6 x 10 8 M Mpc" 3 . 

Using 2dF, iMadgwick fc et al.J (|2002h and lCroton fc et al.J 
(2005) found that 



Mpc~ 



Thus, the mean (M*/Lb) is ~ 4. Given that 80% of the 
baryons in galaxies are in stars (|Read fc Trenthamll2005h . a 
reasonable baryonic mass-to-light ratio is (Mh&r /Lb) ~ 6. 

Table IB1I shows the parameters used along with their 
contribution to the error budget. 



B2 Metals lost at z = in terms of stellar masses 

In terms of stellar masses, Eq. [21] then becomes in terms of 



M G , lost (M.) = 12§ (M*)M* 
ri 



t/eff(M* 



Mr. 



(B3) 



The relationships needed here are the stellar mass- 
metallicity (M*-Z), the y e tt(M+), th e stell ar TFR, and the 
stellar mass function. iTremonti et alj (|2004[ ) constructed the 
M+-Z relation from 53,000 SDSS galaxies and found: 

12 + log % = -1.49 + 1.85 (log M*) - 0.08 (logM*) 2 . (B4) 
H 

They also modeled the effective yield as 

Vo " 21 



log - — = — log 

Vo 



(Vo\< 

\v Iot J 



(B5) 



10 It should be noted that M in the mass-to-light ratio M/Lg 
should be the baryonic mass since the factor yo /y e ff — 1 estimates 
the fraction of ejected metals, as a fraction of the total baryonic 
mass of a galaxy. 



where yo = 0.0104 and Vp = 85 kms We u sed the stellar 
TF relation M*-V rot of lBell fc de Jond (|200ll ): 



AL 



10 9 ' 8 M Q 



( Z V 
V 100km s" 1 / 



(B6) 



The final ingredi ent is the stellar mass function from 

iRead fc Trenthaml (|2005l ). 



B3 Metals lost at z = 2.5 

For our calculation at z = 2.5 in sect ion 13.31 we used the 
same y e g relation from iGarnettl (|2002tl fEq. IBT). Whether 
the 2/ e ff-Vrot relation is the same at those redshifts is 
still an open question. However, one can argue that this 
is likely for the following reasons: (i) the 2 = y a a- 
Vtot relation is made of galaxies of very different star for- 
mation histories (i.e. forming over a wide range of red- 
shifts), (ii) V ro t is a measure of the potential well, which 
is a key factor determining the outcome of th e galactic 
super- wind phe nomenon (e.g. lAguirre et al.l 20011). and (hi) 
both local (e.g. iLehnert fc Heckmanl 1996al : Heckman et alj 
| 2000l: iDahlem et al.lll997f) and high-re d shift starbursts (e.g. 
IPettini et al.l l200ll : fAdelberger et al.l 120031 ) drive strong 
galactic winds. 

We note that t he pa rameterization of y c B-V I0 t relation 
by ITremonti et al.l (|2004l ) (Eq. IB5[I diverge at the low mass 
end, and our result would increased by 25%, a factor that 
depends strongly on the minimum luminosity M m i n used in 
the integration of ■ We use d the parameterization of 

Ve«~Vrot relation of lGarne ttl ( |2002f ) fEq, IBl[l . which does not 



diverges and thus produces results that do not depend on 

the minimum lu minosity M m j n . 

For the LF, Isteidel et ail l|l999f ) constrained the faint 
end slope of the 2 = 3 LF and found it to be much steeper 
a = —1.60 than the local LF. H owever, recent detailed work 
bv lSawicki fc Thompson! (|2006l ) did not confirm such a steep 
LF. They found a = -1.43±0.15 a t z = 3 and a = -1.2±0.2 
at z — 2. Gial longo et al.1 (|2005l ) studied the evolution of 
the LF in a parametric way, and kept the faint end fixed 
to a = —1.3 with redshift. They found that M* ° (2 = 

4400A.AB 

2.2) = -20.93 - 1.48 log(l + 2.2) or M* = -21.67. 

4400A.AB 



ISawicki fc Thompson! (2006) found that M* 700 ab = —20.60, 
which corresponds to M^ i00 AB = —21.60 where the K- 
correction was —1.03 assuming a SED slope of /3 = —1. Both 
studie s found similar <j > *: </>* = 0.0026 (4>* = 0.003) accord - 
ing to iGiallongo et all l|2005l ) l|Sawicki fc Thompson! I2OO6I ). 
We converted the numbers for a h = 0.7 ACDM universe 
when necessary. 

W ith regards to the mass-m e tallicity or luminosit y re- 
lation. iKobulnickv fc Kocl (|2000l ). IPettini et all (|200ll ) and 
iMehlert et al.l (|2002h have shown that 2 ~ 3 galaxies are 
2-4 mag over-luminous for their metallicities, when com- 
pared with the local m etallicity-luminosity r elation. A con- 
clusion also reached bv lShaplev et alj l|2004l ) at 2 ~ 2 on a 
small sample of 8 galaxies. Conversely, at a given luminos- 
ity, galaxies a re a factor 2 to 3 more metal poor. Following 
on this work, lErb et al.l l|2006h used a much larger sample 
of 87 2 ~ 2.2 galaxies and constrained the stellar mass- 
metallicity (M*-Z) relation. They found that 2 ~ 2.2 UV- 
selected galaxies are 0.3dex less metal rich at a given stellar 



The missing metals problem. Ill 17 



mass than the z = sample of iTremonti et al.l (|2004l ) F^l . 
lErb et all (|2006l ) also found that the scatter in the M+-Z 
relation is much less than in the Lb-Z relation, owing to 
the larger scatter in mass-to-light M/Lb ratios. 

The evolution of the TFR to z ~ 2 (10 Gyr 
ago) is unknown, and its evolution to z = 1 is de- 
bated. Whether it evolves or not appears to depend 
mostly on the photometric band. In summary, the B- 
band TFR s eems to be offsetted from the local TFR by 
about 1 mag jMilvang- Jensen et al. 20031 ; Bohm et al. 2004 



Table B2. Parameters used in section [3.31 . 



Nakamura et all 120061 ; iBamford et all l200fj ; IWeiner et all 
20061 ). However. IVogtl (|200lh argued that most of the evolu- 
tion is due to selec tion effects. The if-band TFR evolutio n 
appears negligible l|Conselice et al . 2005; Flores et al . 2006). 
Given the uncertain conclusions, we will first assume an off- 
set of 1 mag for the z ~ 2 TFR, and repeat the calculation 
with no evolution of the TFR. 

The mass-to-light rati o M/Lb is mu c h dif ferent at 
z = 2-3 than at z = 0. iBell fc de Joiigl i|200lf ) showed 
that M/Lb is a strong function of color and is about 1.5 
for the bluest (youngest) galaxies. Similarly to the pre- 
vious z — calculation, we can take the ratio between 
Jb a nd pi, at z ~ 2.2 to estimate the mean (M/Lb) 
ratio. iRudnick et al.J (120031 ) and o ur Fig. [T] showed that 
p*(2 = 2.2) ~ 2 x 10 s M Q Mpc -3 . IRudnick et ail (|2003l ) 
estimated js(rest) at z = 2.2 and found that .7 b ~ 10 s 8 3 . 
Thus, (M /L B ) is about 1 to 2. lErb et all [|2006h showed that 
the M/Lb ratio for UV-selected galaxies at z = 2, ranges 
fr om 0.02 to 1.4. W e will use a stellar mass-to-light ratio of 
l. lErb" et al. I i|2006h also indicated that these galaxies have 
large gas factions, from 20% up to 80%, a result in con- 
trast to local galaxies where the gas fractions are ~ 15% on 
average (see Table [Cl) . Thus, a reasonable baryonic mass- 
to-light ratio for high redshift galaxies is 2. 

Table IB2I shows the error analysis for each of the pa- 
rameters coming into play in our calculation. One sees that 
the dominant uncertainties comes from y e g and the mean 
(M / L) ratio and our result (Eq. l28p is uncertain by a factor 
of ~ 2. 



Parameter 



Value Error Relative Relative 
(1-tj) Error Error to po 

(%) (96) 



2/eff 

LF 0* 
LF Mi, 
LF a 
TF slope 
TF intercept 
TF offset 
LZ slope 
LZ intercept 
LZ offset 
M/L 



0.003 
-21.69 
-1.3 
-7.3 
-20.11 
-1 

-0.185 
5.238 
-0.30 
2 



0.001 
0.25 
0.2 
0.5 
0.25 
0.25 
0.01 
0.2 
0.07 
0.67 



Uncertainty to Eq. [28] 



25 
30 
1.1 
15 
7 
1 
25 
5 
4 
25 
30 



+1.1E6 
-1.2E6 
+7.4E5 
-7.4E5 
+2.9E5 
-2.7E5 
+1.4E6 
-7.3E5 
+1.4E5 
-1.2E5 
+4.0E5 
-3.5E5 
+4.0E5 
-3.5E5 
+1.2E6 
-8.0E5 
+1.3E6 
-8.2E5 
+4.5E5 
-2.8E5 
+1.1E6 
-1.1E6 



+2.9E6 
-2.2E6 



APPENDIX C: METAL BUDGET 

In this Appendix, Table [ClI summarizes the metal budget. 
Some classes of objects are repeated and care should be 
taken in combining different entries. 

This paper has been typeset from a Tp^X/ IATpjX file prepared by 
the author. 



lErb et al. (2006) converted the lTremonti et al. I j2004l) metal- 
licities to the same [NII/Hcr] scale. 
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Table CI. Contributions of various classes of objects. We use Zq = 0.0189, and in our cosmology, the critical density is p c = 1.36 
10 11 h'j M e Mpc" 3 . 



(M Q Mpc" 3 ) 



p/Pc p/p b Reis 
(%) 



Z/Z B p z Pz/Pc 

(M© Mpc" 3 ) 



Z b Pz/Pz, tot Refs ' 
(Z ) (%) 



Baryons (p b ) 



5.98 x 10 9 0.044 



100 1,2 



z = 



All Metals 










7.75 x 10 7 


5.70 x 10~ 4 


0.6854 




1 


MS stars 










4.35 x 10 6 


3.20 x 10~ 5 


0.0385 




1 


WDs(C+0) 










4.90 x 10 7 


3.60 x 10~ 4 


0.4329 




4 


Metals produced by Type II SNs 








2.30 x 10 7 


1.69 x 10~ 4 


0.2034 


100.0 




Stars z > 


6.43 X 10* 


0.004728 10.75 f SFR 


1/42 b 


2.13 x 10 7 


1.56 x 10~ 4 








Stars 


5.63 x 10 s 


0.004140 


9.41 1,2,5 (Salp.) 


1 


1.48 x 10 7 


1.09 x 10" 4 


0.1307 


64.3 




BH+NS 










1.60 x 10 7 


1.18 x 10~ 4 


0.1419 


69.9 


4 


Hi 


5.71 x 10 7 


0.000420 


0.95 8, 9 














H 2 


2.17 x 10 7 


0.000160 


0.36 10 














Cold gas (+He i) 


1.06 x 10 s 


0.000782 


1.78 (8+10) x 1.35 


0.84 


1.69 x 10 6 


1.24 x 10~ 5 


0.0149 


7.3 




Hot gas (cl.) 


3.57 x 10* 


0.002625 


5.97 3 
















2.45 x 10 s 


0.001801 


4.09 4 


0.3333 


1.54 x 10 6 


1.13 x 10~ s 


0.0136 


6.7 


10 


Hot gas(gr.) 


6.80 x 10 8 


0.005000 11.36 3 


0.3333 


4.28 x 10 6 


3.15 x 10 -5 


0.0379 


18.6 




WHIM IGM OVI 


> 3.67 x 10 s 


0.002700 


6.14 11 


0.1? 


6.94 x 10 5 


5.10 x 10~ 6 


0.0061 


3.0 




WHIM IGM BLyAb < 7.89 x 10 8 


0.005800 13.18 11,18 


0.1? 


1.49 x 10 6 


1.10 x 10~ 5 


0.0132 


6.5 




2 = 2.5 


Stars 2 < 2 < 4 


1.21 x 10 8 


0.000890 


2.02 J SFR 


1.260 


4.00 x 10 6 


2.94 x 10~ 5 


0.0354 


100.0 




BX 


4.02 x 10 7 


0.000185 


0.67 


0.5 


3.80 x 10 5 


2.79 x 10~ 6 


0.0034 


9.5 




BX+K20 


1.98 x 10 7 


0.000146 


0.33 


0.8 


3.00 x 10 5 


2.21 x 10- 6 


0.0027 


7.5 




DRG 


1.06 x 10 7 


0.000078 


0.18 


1 


2.00 x 10 5 


1.47 x 10" 6 


0.0018 


5.0 




SMGs > 3mJy 


2.73 x 10 7 


0.000201 


0.46 


<1.86 


<3.63 x 10 5 


<2.67 x 10~ 6 


<0.0032 


<9.1 




DLA (N > 20.3) 


1.50 x 10 8 


0.001103 


2.51 13 


0.07 


1.98 x 10 5 


1.46 x 10~ 6 


0.0018 


5.0 




30% missed 


5.00 x 10 7 


0.000368 


0.84 14 


0.7 


6.62 x 10 5 


4.86 x 10~ 6 


0.0058 


16.5 




LLS 19-20.3 


3.60 x 10 7 


0.000260 


0.60 19 


0.1 


< 6.80 x 10 5 


< 5.00 x 10~ 6 


<0.0060 


<17.0 


19 


LLS 19-20.3 


2.72 x 10 7 


0.000200 


0.45 19 


0.15 


> 7.71 x 10 4 


> 5.67 x 10~ 7 


>0.0007 


>2.0 


19 


LLS 17-19 




















Forest 


5.44 x 10 9 


0.040000 90.91 3 


0.003 


3.08 x 10 5 


2.27 x 10~ 6 


0.0027 


7.7 




Forest 


6.14 x 10 9 


0.045130 102.6 


<0.005 < 5.80 x 10 5 


< 4.26 x 10~ 6 


<0.0051 


<14.5 


19 



Refer ences : (1) iKochanek et alj feOOll), f2llBellet al.l J 2003h. (31 IFukueita et al.l dl998ll, (4) iFukueita fc Peebles! (l2004fl . (5)[Coleet_ 

et al.l ||20031. ( 9 1 [Rosenberg fc Schneider! j2003h . (101 iKeres et 



l|200lh . (6llBell et alj ll2003h. (71lRudnick et a, 
"bach etjdj J2004h. (12 1 



000. 



1 



11 



.1 <l2003h. (8) Zwaan 
d2003h. 



Dunne et al 



Simcoe et alj l|2004l l. (16l lSchavel l|2001bl l. (171 



13 



Pettini et al 



iPeroux et alj d2003f). 
(|2003h . 



114) 



(18l lRichter et al 



Vladilo fc Perouxl <2005n . (15) 



2006 



(19) This paper. 

'Averaged yield for type II SN with m > 10 Mq llMadau et al.|[T996ll . taking into account a recycled fraction fo R = 0.28. 



